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A. BACKGROUND 


As combat systems aboard naval surfac2 ships continue to 
grow in importance, size, and weight, the feasibility of 
their installaticn aboard lightweizyht vessels can only be 
realized by reducing other major shioboard weight require- 
ments. A significant achiavemen* in reducing main propulsion 
weight was reached with the advent o£ gas-turbine propulsion 
plants, of which the majority of our new combatants will be 
powered. However, for the remaining capital ships and subma- 
Gem=s £O be constzructed using 2ither conventional or nuclear 
st2am systems, power plant weight reduction can only be 
realized by increasing the effectivenass of tha individual 
plant compcnents and thereby reduciag their sizes. Plans call 
fee couture installaticn of Ramkine-sycle waste-heat recovery 

u 


? 


opulsion, so 
Significant weight reduction could also be realized there 


systems even on those ships with gas-turbine pr 
through the use of more effective systems. 

Mice cOmeche Nature Sf ehe condensing process, repoyes 
greatest thermal resistance to condenser-tube heat transfer 
occurs on the tuhe side, anda thorough review of tube-sidé 
enhancement is presented by Bergl2s and Jensen [Ref 1}. 
Internal enhancement of the tube without regard to the 
Shell-side resistance problem, how2vear, would be a wasted 
effort if improvements were made to tha point where the 
external (shell-side) resistance became the controlling 
Baes+Or in the transfer of heat. Recent investigations into 
the shell-side condeéensing process and techniques to improve 
this process are thoroughly reviewed by Marto [Ref. 2}. 


Outside enhaacement techniques include the use of low 
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Bepegdsal ©£2nS, toped tubes, fluted tubes, and 


fu 
‘CO 
of 
fee 
j s 
iD 
Qu 


Saas angqs tO promote drepwise condensation. 

Ongoing research is underway at the Naval Postgraduate 
Schocl to further study tha shell-side enhancement of marine 
condensers. An endurance éoparatus iS in opération to 
examine the dropwise-promoting effectiveness of various 
metallic and polymer coatings over prolonged periods of 
UM e « Reilly {RBRef. 3] conducted i study of imprevaments 
using various spirally fluted tubes under flimwise condensa- 
mon, and «he effects of tube bunila inundation have also 
been studied by Kanakis {Ref. 4]. 

A test apparatus has been constructed by Krohn [Ref 5] 
to systematicaliy study steam condénsazion on a single hori- 
zontal tube. The apparatus is modeled after a similar design 
Constructed by Rose (Ref. 6], but can also operate under the 


Pegi-Vvacuum conditiors found in Marine condensers and 
without the presence of noncondensable gases. Ehis Les< 
apparatus was instrumented and tasted by Graber [ Ref. 7], 
who determined the Sieder-Tazte coefficient for the test zube 


memgech to b2 0.029. 

Im @2Xperiments using low intégrai fins at Queen Mary 
Mollece of London, Rose, et. Sige tj i. 68.) OCbtargaed an 
wermyum pitch of 2.0 mm (0.08 in.) for a given fin size of 1 
mome.04 in) high and 0.5 mm (0.02 in WeCeti-2O= Tests 
conducted at 9ne atmosphere. Tecwecnes Dp kech, hea reported 
Vapor-side heat-transfer enhancements of 400% and 309% for 
heaz fluxes of 0.3 MWyme? (9.51x10% Btu/hr-ft2) and 0.8 MWyme 
(2.54x105 ctu/hr-it 2), respectively, and for a vapor 
velocizy of 0.7 m/s (2.3 ft/sec). Further testing is needed 
*o compare these results with Similar fin geometries under 
vacuum conditions, while varying fin dimensions (height and 


Sagem) in addition te the fin pitch. 


yZ 


B. OBJECTIVE 


vt 
54 
MD 
i 


The main objectives of this research effort were, 


1) ensure a vacuum-tight apparatus so that data could be 
taken at both atmospheric and vacuum conditions with no 
detrimental effects due to the presence of nonconden- 
Sable gases; 

2) cake data for an instrumented smooth tube to verify 
*he Siader-fate coefficient obtaiaed by Graber; 

3) take data for a smooth tube to check the reliabilty 
of the apparatus and the data-reduction procedures used; 
and 

4) take data for six externally-enhanced tubes of 
Maes ous £in pitches to obtain the relative optimum pitch 


mor a fixed fin gecmecry. 
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A. SYSTEM OVERVIEW 


The apparatus used for this research was essentially the 
same used in references 5 and 7, with several noted modifi- 
Sections. Rh SChematac Sketch Of the system is shown in 
Figure 2.1. Steam was ganerated in a 304.8-mm (12-in) dian- 
eter Pyrex glass boiler by ten 4000-watt, 480-volt Watlow 
immersion heaters. Passing througa a 304. 8-ma (ie —2n) tS 
152.4-mm (6-in) reducing section, the steam travelled upward 
mneOough a Pyrex secticn 2.44 m (8.) £t) in length, around a 
180-degree bend, and back down a straightening section 1.52 
m(5.0 ft) in length before entering the stainless-stesi 
test section. The condenser tube to be tested was mounted 
horizontally in the test section behind a viewport to permix 
visual observation of the cond2nsinjg process. 

Steam that did not condense on the test tube passed into 
@a stainless steel auxiliary condenser, and all condensat:; 
waS returned via gravity to the boiler. Thea auxiliary 
condenser Was NeCOUSE=RUGCtTed Of tWO 9.5-mm (3/8-in) water- 
eooled ccpper lines helically coiled to a height of 457 ma 
(18 in). 

Cooling weter for the test cube was provided by two 
centrifugal pumps connected in series. Th water could be 
@meOttied from zero flow to 0.69 I1/s (11 gpmjy. The maxinun 
water velocity which could be obtained through the tube was 
5-48 m/s (18 ft/sec). A continuous suppiy of tap water was 
used for cooling the auxiliary condenser. Theot: lerg the 
flow of tap water through the condinser was the means used 
to vary the internal pressure of the test apparatus. The 


water flow through both the test tube and the auxiliary 
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need): 
valves and measured by rotameters with full-scal¢ ranges o 


0209 l/s (11 gpm). 
Miata €jeGeon  DECViaGed  £Or removal of noacondensable 


condenser was regulated by 19.1 mm (3/4 in) dianeter 


b 


rh 


gases from ‘the auxiliary condenser through a 12.7-mn 
(1/2-in) line. The source for th2 air 3jector was 1.1-MPa 
(169-psig) house-air supply. 

Be SYSTEM INSTRUMENTATION 


The input voltage through the h3atears was varied through 


@ panel-mounted potentiometer. 440-VAC line voltag2 was 
mecauced by a factor of 100 when fed into a LEfeEr ential 
input precession voltége attenuator. The stepped-dewn 


voltage passed through a True-Root-M¢ean-Square converter 
stage on which the integrated period was reduced tc anourt 1 
ms. The output of the TRMS converter was then bufted and 
compared to a refarence voltage from the potentiometer. The 
Semparator output was fed to the control input of a Halnar 
silicen-controlled rectifier power supply which appvlied «he 
actual voltage to the heaters. The ITRMS cconverter output was 
also paralleled to a filter an than input to the data 
Seguisition system. This input was proportional to «he power 
supply output. A diagram of *he system is shown in Figure 
ea Le 

The internal pressure of the system was measured manu- 
ally by a U-tube, mércury-in-glass mnanomet27 graduated in 
millimeters. Unavoidably, steam could condens2 in the manon- 
Seen. therefore, the varying height of the water column in 
the mancemeter needed to be accounted for when measuring the 
system pressure. 

Temperatures +threughout the systsm were aneasured by 
Sepper-constantan «hermocouples: six for ‘the wall of 42 


speciaily-ccenstructed test tubs, «wo for the steam, and one 
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Pach) foe the cooling-water inlet, condensate return, and 
m 


ambient. The calibration procedure for these tharmoccuples 
is described in Appendix A. The an MiacwrIseeenrOughn 216 
test tube was measured by a Hewlett-Packard (HP) 2804A 


quartz thermometer. 
All temperature measuraments were fed diractly into th: 


data-acguisition system as described below. 


C. DATA ACQUISITION 


mete SG97A Data ACgQuLSition/Control Unit 


W u 
monitor system temperaturas. This was interfaced with an HP 
ae e 


9826A computer which served as a controlling unit through an 
interactive G@dua-Eeduction — progra i ard user keyboard 
prompts. Raw data gathered by the data-acquisizion svsten 
were stored on computer disks for later reduction and evalu- 


aeLon. 


D. SYSTEM MODIFICATIONS 
1. Boiler 


The fiberglass insuiation was removed from «he 
boiler to allow the operator to nore easily monitor he 
water level. Although a closed-system design was used, it 
WaS Still possible fcr sta@am to estap2 via the air ejector 
or through the reliefé valve [Fig. 2.1]. Calculations showed 
that the additional heat loss dua to the removal of «the 
insulation was minimal { App. 8B], and =hesauthor Sele this 
loss was much more acceptable than risking damage tc the 
immersion heaters through a low-water casualty in the 
boiler. 
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2. Condensate Piping 


As originally designed, draining the system required 
breaking down the condensate piping. While this was not a 
daily cccurrence, the procedur2 was inconvenient and the 
existed also the possibility of losing the vacuum integrity 
of the system each time it was dons. 

To avoid these problems, the existing fill-line 
valves were rearranged as shown in Figur2 2.1. This 
arrangement also added two additional features to the 


system: 


1) the fill/drain valve could be spened during operation 
to drain any heavy particulate matter from th? systen - 
Similar to the "bottom blow't procedure used on Naval 
boilers; and 

2) after extended periods of inactivity while opened tc 
the atmosphere, the entire system could be given a thor- 
Sugh steam-cleaning by following the procedures outlined 


ie ADPpeENndix C. 


3. Vent Valve 


Piemewodditcatson Of the sondéensats return piping 
necessitated the addition of avant valve for use when 
filling or draining the system. A &Y.3-mm (0.17-in) needles 
valve was installed on the 101.6-um (4-in) flange of the 
test section. This valve would also serve as «he tap for the 
proposed sampling of noncondensabi2 gas concentrations in 
the systen. 


u. 
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The original system design used a 6.4-mm (1/4-in) 
Stainless steal tube angled down t> the mercury manometer. 
During this thesis, the manometer was raised +9 eye level to 


facilitate easier and quicker reaiing. Replacement of the 
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stainless-steel line by 2 12.7-mm (1/2-in) ccpper <ub2 
reduced the possibilty of error caused by water slugs 
building up in the smaller diameter tube. The mere workable 
copper was chosen over stainless str2eal to reduce the stiff- 
ness of the connecting line. This was necessary to 
eliminate leakage in this part of the apparatus as explained 
maesoection II.£. 


As an alternative +o the manometer, a cCelesco 
Strain-gage pressure transducer was installed on the tesr 
section flange next to the vent valve. The calibration line 
Mor the transducer is shown in Figure 2.3. The author felt, 
however, that the reliability of this measurement would not 
be high enough until a second, more accurate transducer was 
installed. Once inccrporated into the system, theugh, “hese 
transducers would provide automatic input to the data acqui- 
Sition system, Eliminating the requirement to manually enter 


the manometer reading into the data-réeduction progran. 


6. Relief Valve 


AS originally designed, 2 6895-Pa gage pressu 
(1.0-psig) relief valve was installed beneath a 1.0-m 
(39.4-in) Length of 12.7-mm (1/2-in) stainless-steel pi 
Steam which condensed and became trapped in this piping 
would open the valve at a water-ccoluma height of only 0.70 n 
(27.6 in). Once opened, a back-pr2ssure of 0.14 MPa gags 
(20 psSidq) was required to réseat the valve - something unob- 
tainable even with an absolute vacuum on the inlet sids of 
the valve. To avoid this problen, the valve was raised in 
the line to a point only 76 mm (3.0 in) below the outlet 
from the auxiliary ccndenser section and was replaced by 
another 6895 Pa gage pressure (1.0-psig) relief valve which 
reseated at only 0.04-MPa gage (6-psig) back-pressure. 
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The steam temperature prope was located directly 
above the test tube and shared th2 same inlet to the te 
section as th2 manometer line. When drawing a high vacuu 


n 
this arrangement allcewed water whith had collected in the 


(tt 


Manometer to be drawn back into <he test section where i 
would flow down the probe and onto th? tube. To prevent this 
water and any contaminants picked up in the manometer from 
being deposited on the tube, the probe was bent sce that it 


mee CLiset from the center or the t2ast tube. 


8. Cooling Water Sy 


GP ef 42 as eS ae 


wn 
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A second centrifugal pump was added in series to the 
one already installed. This pump boosted the maximum cooling 
Beeert velocity through che test tube from 3.96 m/s (13 
fr/fsec) to 5.48 ms (18 ft/sec). 

Meme Some (5) £Om) Aalr=cOnditioOning unit used wit! 
the cooling water system would enérgize at a water “*2mpera- 
ture of 17 °C (62.6 °F) and sacure when the <anperature was 
meauced to 13.4 °C (56.1 °F). Therefore, for a given stéan 
Bemoperatura of 50 °C (122 °F) azemmeene = ube, he log- 
mean~temperature difference would vary by as much as 11%. 
But the measured temperature rise of the water through the 
mmoe Showed very little change. fo avoid this transient 
problem, the alr-conditioning unit was not used, and instead 
fresh tap water was continuously f2d to the sump while an 
equal amount of water was being drained, maintaining a 
constant sump level. This methoi provided a constant- 
temperature supply of cooling water +9 the inlet side of the 
test tube. 

irovoirwonn= poOlinzed Out that, while the duty cycle 
Memeene air-conditioning unit was 2 function of «he ther- 


mostat used, a more sensitive thermostat would require the 





use of a more expensive cooling device than a commercial air 
conditioner. An alternative solution would be to install a 
much larger sump in the system. 


9. Thezmopile 


Since the temperature ris2® of th? coslin water 
through the tube was a critical measurement in the 2xperi- 
ment, a 1l0-junction thermopile was added to neéasure this 
temperature rise in addition to the quartz thermometer. As 
will be explained in Section V.A.3, howaver, problems arose 
With the thermopile, and it couli not be accurately used 


M@irang this thesis. 


E. VACUUM INTEGRITY 


One of the objectives of this thesis was to ensurée a 


h 

Vvacuum-tight test apparatus to eliminate the presence of any 
noncondensable gases and their detrimental effects on the 
condensing process. A standard of no more than a 5.0-mnm 
mercury (0.10-psi) loss over a 24-hr period was considered 
to be an acceptable tolerance, but obtaining a leak rate 
within this tolerance proved to be tne nost time-consuming 
effort during the research. 

Due to th2 construction of the apparatus and nature of 
the experiment, most leak-detection methods could net bé 
used. A sealing substance could not be used without risking 
Contamination of the interior of the apparatus which might 
prohibit filmwise condensation on tne test tube. 

Initially the system was pressurized and the standard 
soap-soluticn test was used to locats leaks. Once pressur- 
ized, a liquid-soap sclution was applied to each external 
joint or fitting where a leak could be present. The higher 
pressure air inside the apparatus would escape through any 


leaks and produce bubbles on the applied soap film. 
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However, the maximum pressure which the Pyrex glass members 
Gould tolerate was Only 0.074 MPa gage (10.7 psig), and for 
safety reasons the authcer chose not to pressurize the svstem 
to more than 0.034 MFa gage (5.0 psig) - a pressure differ- 
ence between the atmcsphere and th2 apparatus of only 0.040 
MPa (5.8 psi). The numerous 2xtecnal valves and fittings 
prohibited the use of an evacuated hood to achieve a greater 
rressure difference. 

A similar test could not be used to locates any vacuun 
leaks which were not present when th? system was under a 
positiva pressure, as the apparatus was not larg¢ enough to 
permit «he apolication and observation of a@ soap solution on 
the interior. 

In another attempt to locate the leaks, a National 
Research Corporation (NRC) 101.6-mm (4-in) vacuum pumping 
system was connected to *ha apparatus. This pumping system 
included a Welch model 1376M mechanical pump and a model 
NHS-4 diffusion pump. An NRC model 521 thermocouple gage was 
aiso connected to the test apparatus and the entire appa- 
ratus was evacuated to 0.21 torr (4%.1x1073 psia). Acetone 
was sprayed around all flanges, fittings, and joints. A leak 
around any of these should have produced a rapid rise in the 
thermocouple reading, but this method also proved ineffec- 
eave » proktably due to the large size of the apparatus 
resulting in too great a mean free path for the acetone 
molecules to travel from the leak t5 the thermocouple. 

The next alternative was +o braak the system apart into 
three main sections: the glass boiler and steam piping, the 
stainless-steel test section and auxiliary condenser, and 
the condensate return piping. The glass section was blank- 
flanged and evacuated to an absolute pressure sf 0.033 torr 
(6.4x10 -* psia). The rate-of-riseé measurement for this 
section showed a loss of only 0.48 mm of mercury (0.01 psi) 
over a 30-hr period. It was, thereforea, concluded that any 


leak in the assembled apparatus was not from this section. 
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Once removed, tiem eects ece On and the auxiliary 
condenser were blank-flanged, pressurized to 0.10 MPa gags 
(15 psig), and immersed in a large plexiglas tank filled 
with water. This test easily revaaled anumber of smali 
leaks abcut the inlet sid= of the test tube and also around 
the plug which connected the condensate return piping to the 
Pease Of the auxiliary condenser. Rsplacing an O-ring in the 
meuape fitting and siiver brazing the plug into place elinmi- 
nated these leaks. 

The same immersion t2est revealed small leaks in the 
joints of the condensate return piping. These leaks were 
eliminated by replacing all stainless-steel farrules in the 
Swagélick fittings with teflon ferrules. 

Once reassembled, considerable leakage was still indi- 
cated by a substantial overnight rise in th? manometer 
level. The author felt confident that this leak was not in 
the main assembly of the apparatus, but was instead in the 
Manometer assembly itself sinc2 it could not be pressurized 
for testing. As mentioned in section II.C.4, the stainless- 
meek line leading to the manomate> was at this time 
replaced by a 12.7-mmn (1/2-in) soft copper tube. This tube 
eliminated the need for two 90-degree elbows and three 
lengths of stainless stesl tubing in the line - an assemb 
which proved too rigid to allow even the slightest misalig 
ment into the manometer. 

Upcn completion £ the installation of this assembly, 
the system was evacuated to an absolute pressure of 92.5 mon 
Hg and cver a 24-hr period the mercury l¢ével rose to only 
94.0 mm. This leak rate was well within the acceptable 
tolerance. 
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Le TUBES TESTED 
1. Instrumented Tube 


An instrumented tube was fabricated from a2 thick- 
faeeisd Cepper tibe with an inner diameter of 12.70 Bm (1/2 
in) and an outer diameter of 19.05 mm (3/4 in). The tube 
was cut into three sections ints which six holes were 
drilled axially along th2 walls at equal spacings 60 apart. 
These passages were fitted with 0.094-mm (3/32-in) OD capil- 
Maries { Fig. 2.4] which were silver-soldered ints place, and 
the three sections of tube were than soidéered back into one 
piece. Thermocouples w2re fitt3d into the capillary 
sections to measure an average wall temperature. 

By knowing the average wall temp¢srature, the Nussel* 


humber for the inside could be computed. By computing the 


iD 


gradien+ of the Nusselt number agains<c the Leder-Tat 
parameter, the inside coefficient could be obtained as the 
inverse cf the gradient. Baci@e 2.5 SHOowSs a2 photograph of 


the instrumented tube with the installed thermocouples. 
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A smooth tube with no wall thermocoupl?s was also 
tested +o obtain the inside heat-transféer coefficient 
through the us¢ of a jHodtatederNttson  Ploce yf Ret. 98]. 
Determination of the inside heat-transfer coefficient was 
Sercical to the experiment, as it was used to obtain the 
outside condensing coefficient for the smooth tube and all 


of the enhanced tubes. 


3. Finned Tubes 


TOME iemene W0rime OD jeCr~a vs OL this thesis, 3 
Ser-es cf six finned tub?s was also tested [ Fig. DOr 
These tubes had the same overall dimensions as those above, 
but were enhanced with radial fins 1 mm (0.04% in) high and 1 
momo, 04 injecnicks Each tube had 2 ditferen* fin pitch and 
was tested to determine a relative cptimum pitch. Fin 


Pemeelies tested were 1.5, 2.0, 2.5, 3.0, 5.0, and 10.0 mn. 


G. SYSTEM OPERATION 


The tube to be tested was cleaned in a warm soluticn of 
Sparkleen and then rinsed with tap water, which produced a 
contaminant-free, wetted surfacs. The tube was then 
eiscadled in ths test section, cars being taken not to touch 
Be cOncaminete the ccndensing surface. 

The systen was) UEOUGhe. tO  Opeeating pressure by 
following the procedures of Appendix D, and data collection 
began when steady-state Conar= 200s were achieved. 
Steady-state conditicns ware determined by observing the 
steam temperature measured by the respective thermocouplss. 
When their output voltage on the HP 3497A reached a ccnstent 
Merue With fluctuations of only on? or two microvolts, it 
was assumed that steady-state conditions existed in the test 


apparatus. 
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Data sets were taken starting with a test-tube cooling 
Weer tlow fate of 90% (Which corresponded to a water 
Memocicy Cf 4.95 m/s or 16.2 £t/ses) , ranging downward in 
decrements of 10% through a minimum flow of 20% (1.16 m/s, 
3.8 ft/sec), and then upward from 25% to 85% in increments 
of 10%. After adjusting the flow rats, the tempezraturs rise 
through the tub2 was monitored by observing the digital 
output of ths quartz thermometer. When this rise became 
Gonstant, another set of data could be taken. Duriag data 


decrease in 


Dp 


Bims, a siight rise in préssure accompanied th 
meme cooling water flow through the tube, a resul* of the 
meauced heat flux. Similarly, Hineao@nd sane Llow sh=ough 
the tube caused a slight decreas2 in the system pressure. 
This variance could be anticipated, and Sinc? data were to 
be collected at a constant pressur2, it was easily cempen- 
mmremetor bY throttling the flow of cooling water through 
the auxiliary condenser one or two percentage points on the 
rotameéeter. 

Something which cculd not be anticipated, however, were 
Sudden rluctuations in the tapwater pressure to the auxil- 
lary condenser which caused pressure changes cf saveral 
Millimeters o£ mercury inthe systan. Loma voe de oc his 
problem, the flow through the auxiliary condenser had +o be 
continuously monitored, unless data were being taken late at 
night when there was no demand on the laboratory building's 
Water supply. The test tube could be easiiy monitored 
Mmeough the viewport for confirmation of filmvise condi- 
fons. If there was any sizeable change tc dropwise 
condensation on the tube, the data set was discarded and the 
procedure was repeated. 
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III. FILMWISE CONDENSATION 


—_ ae aw ee ee ee ee 
e 


A. THE DROPWISE PROBLEM 


It was essential during the course of this thesis to 
collect data under filmwise conditions. Numerous problems 
wer2 encountered by Graber [{Ref. 7] in avoiding the *ran- 
SresoOn tO dropwise condensation during operation, and his 
proposed solution was a vacuum-tight test apparatus. Since 
the tube surface would wet completely after installation, 
contaminants leaking into the system were possibly adhering 
to the surface and promoting the dropwise condensation. 

However, even after »btaining a vacuum-tight apparetus, 
the author was still unable +o maintain good filmwise 
condensation for more than two hours on a smooth tube. While 
this was enough time to collect a complete set of data for 
this tube, filmwise ccndsnsation last2d seventesi minutes at 
the most on any of the finned tubes, the average time being 
less than tén minutes. This was vraiictable - the corners of 
the fin/surface interfaces provijied a better trap for 
contaminants and were harder to cl2an - but unacceptabls. 

The use of the stsam-cleaning method described in 
Appendix C would thorcughly clean the tube so that complete 
filmwise condensation was re-established, but dropwise 


condensation would again become prevalent within minutes. 


Be. SOLUTION 


Having eliminated the possibilities of iastalling a 
dirty tube or contamination due to leakage, the only reason 
for the dropwise problem had to be coming from outgassing of 
the nylon holders for the test tube. The outgassing rate for 


nylon was found to be almost two oriers of magnitude greater 


fy 





(D 


than +he rate for Teflon [{Ref. 9], Teflon being the san 
order as stainless steel. By punping down the apparatus 
while it was hot, nyion molecules were being outgassed into 
the test section and were immediately being déposited onto 
the surface of the cooler test tub2?, resulting in inadver- 
tant "sputcering" of the tube with a nylon coating and 
subsequent dropwise condensation. To eliminate the nylon 
mmctcerface with the interior of the test section, special 
stainless steel caps were manufactured to slip over thea 
nylon holders. 

Teflon bushings were fitted within the caps to insulate 
them from the the tube. Teflon had both alow thermal 
Gonductivity and a relatively low sutgassing rate. Figure 
3.1 shows a detailed sketch of the installation. 

This configuraticn appeared to solve the problem of 
dropwise condensation. Although an actual endurance tes* was 
mou Conducted, the system was in operation intermittantly 
for over fifteen hours with ¢he smooth tube and over four 
hours with each fanned tube with n> breakup of the filmwise 
condensation. 

Once installed, this arrangemant also eliminated the 
need for the tube-cleaning procedure recommended by Graber 
(Ref. Vj, who felt that a strong cl2aning solution of sodium 
hydroxide and ethancl was need2i to decontaminazte the 
pur cace. Only a warm solution of Sparklieen was used 


throughout the data-ccliection stage of this thesis. 
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C. VAPOR VELOCITY LIMITATIONS 


To check the outside heat-transfer coefficient data with 
the Nusselt prediction, and also t9 obtain @n accurate 
inside coefficient wit the smecdidtied Wiison Pict, vapor 
velocities approaching zero were preferred. Due to wzhe 
design of ‘the system, however, the pressure drop to the 
auxiliary condenser required vapor flow past the test tube. 
This being the case, attempts were made to minimize this 
flow velocity by cutting down the power to the boiler and 
mBeGcottiing back on the ccoling water supply tO the 
condenser. AS vapor velocity was dscreased, however, drop- 
e 


*seast tube. Under 


rea] 
Er 
ae 


wise condensation again took place o 
axmospheric pressure, this occurrence was at vavor veloci- 
tias of about 0.5 m/s (1.6 ft/sec), and under vacuun 
MeesationS iit occurred at apout 0.9 m/s (3.0 tt/sSéec). 
@Meparen=zly there still axisted a sizeable rate of outgassing 
within the test section, those gases collecting about the 
bes. tube and interfering With the filmwise condensation 
process. This suspicion was confirmed when an increase in 
mae Vvapor velocity eliminated the formation of dzecps on the 


tubs. 
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IV. DATA REDUCTION 
The data collected and stored on the computer disks 
reduced uSing the programs WILSON, SIEDER, and ODRP. The 
programs were amenable to changes, which allowed the author 
to analyze and compare results while varying parameters 
within «he programs. Stepwise reduction procedures for each 
program are listed below and the program listings are found 


moe ApDrend:x E. 


Ae. PROGRAM SIEDER 


1. Compute the average cooling water temverature. 

2. Compute the average wall t3anperature. 

3. Compute the cooling water velocity. 

4. Compute the mass-flow rate of the cooling water. 
5. Compute the heat «ransferrai t. the cooling water. 
6. Compute the average inside wall temperature. 

7. Compute the log-mean-ctemperature difference. 

8. Compute the inside heat-transfar coefficient. 

9. Compute the Nusselt number. 
10. Compute the Sieder-Tat2 parameter. 


V1. Compute the inside coefficien:«. 


Be. PROGRAM WILSON 


1. Assume a value for the Sieédar-Tate coefficient. 

2- Compuze the Reynolds and Prandtl numbers for flow 
through the tube. 

3. Compute the log-mean-temperature difference, heat 
flux, and overall heat-transfer coefficient for the «ube. 


4. Assume an outer tube surfac2 temperatures. 
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5. Compute the outside condensing coefficient tusing 
properties é¢valuated at the film tenperature. 

6. Compute the outer surface temperature and iterate 
Beeps 5 end 6 if not within 1%. 

7. Assume a viscosity-correction Faceor “aEcE the 
Sieder-Tate equation cf 1.0. 

8. Compute the inside heat-t*transfer coefficient. 


9. Compute the inner surface temperature. 


fs 
wz) 
ct 
) 
rj 
§») 
ct 
wD 


ieee COMpPUCe che Viscosity correction factor an 
steps 8 threugh 10 if not within 14. 
i. Congute the Sieder-Tate coefficient and iterats 


Beeps 2 through 11 if not within 0.3%. 


Ce. PROGRAM DKP 


1. Compute the average ccoling water temperatures. 
2. Compute the cooling water velocity. 

BeeCOmMpuce the mass-flow rate of the cooling water. 
G. Compute the heat transferred to the water. 


- Compute the leg-mean-temperature difference. 


Compute the overall heat-transfer coefficient. 


Compute the wall resistanc2 of the tube. 


Compute the Reynolds number of the cooling water. 


inside heat-transfer coefficient. 


ct 
ty’ 
(wD 


Compute 


oOo wo ON HD WN 
° 


Compute the condensing heat-transfer coefficient. 
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¥. RESULTS AND DISCUSSION 

Numerous data runs were made using th procedure 
described in Section II.G. Time constraints, nOWéEVELr, 
limited the number cf repeat runs that could be made for 
this thesis. Primary concern focused on establishing 2 reli- 
able, repeatable Sieder-Tate coefficient. Data were taken 
for all of the tubes at both atmospheric and vacuum (88mm 
Hg, 1.7 psia) conditicns. Complete filmwise condensation was 
Maintained for all data runs, and the mass concentration of 
noncondensable gases was held betwe2an 0% and -1% during all 
testing. The negative value was indicative of slight super- 


heat in the system or an inaccurate manometer reading. 


Ae INSIDE HEAT TRANSFER COBFFICIENT 


1. iInstrumented Tub 


iD 
12 6) 


Resulzs 





Figure 5.1 shews the variation of the Nusselt number 
as a function of the Sieder-Tate parameter for the instru- 
mented tube run at atmospheric pressure (run SDA7). This 
method yielded a Sieder-Tate coefficient of 0.035 on two 
separate data runs (SDA7 and SDA8). The same method under 
vacuum conditions yielded a coefficient of 0.037, which was 
also repeated (runs SDA5 and SDA6). The temperature distri- 
bution around the tube wall was symmetrical about the 
vertical plan passing downward through the centerline of 
the tube, and showed as much as a 158 C temperature drop 
meom the top 9f the tube t9 the botton. 


2- Smooth Tube Results 


PeCiEom a,c snows the modified sWalson Plot for 


smooth-cube data ccllected at atmospheric pressure (run 
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ODt2) . Liasmetlea nod produced a Sisder-Tate c 
0.033 with an intercept of 1.16. A Similar pl 
conditions (run 0D10) produced a cozsfficient of 
mererespet Of 1-27. The Nusselt theory predicts an intercept 
mao 5s DUS, aS yeas Noted in Section III.c., a finite vapor 
velocity past the tubé was unavoidable in this resear 
Paice “oss rsnows esa plot Of the Fujii { Ref. 10] 
Berrelation £26 the smooth-tub2 runs. This expe 
Beerelacion is to correct the Nusselt number for = 
of vapor shear on the test tube ani is a plot of the equa- 


On : 
NuRe-!/2 = 0.96F1/5 


The results obtained were Slightly higher chan those 


Peed2cted using “his correlation. 
3. Discrepancies 


tion ore dinamo leden-Tate Squetior [Refs %1) Fer 
meely developed *urbulent flow in 4 tube with an L/D ratio 
Meeater thar 60 has aileading coefficient of 0.027, so a 
Meooer value for the shorter cube (L/D = 18) used for this 
research stands to reason. The data ror the tibes apvears 
Vos, ewes o, 0.0 36) 


man0.037 = shoWing uc to a 6% scatter from the mean. Both 


© rt 


at first inconsistent, with results of 


tubes, however, showed a larger coafficient when the tast 
seczicn pressure was reduced, which reduced the heat flux as 
Beil. Dropping the pressure from atmospheric toca high 
vacuum (88 an Hg, 1.7 psia) reduced the heat flux bya 
Beer Or Of aimost three [Fig. 5.4}. Seeeca enOther Way, 2 
Mowe neat flux produced coefficients of 0.036 and 0.037, 
Mabe @ higher heat flux produced 0.033 and 0.035. 

The author doubts the reliability of the inside co- 


Ber2cient oktained under vactium conditions for two reasons: 


44 





-eqep eqny-qioous roxy [OL “yea ])] UOTRPeTaTIOS Tring €°G einbty 
J 


OOT ul 1 
eS ss oe eee al eateries Te a 


Z1IGO0 UunNY oO nae an 


Otgdo0 uny + \ 


S 1G) Says 


TIf M3 


ap 


42 


"sunT UWNNIPA pre ITAsydsowze ACJ SoaxnTJ eey satzeTSH n°G aernhty 
C2-WyM> b 


OO0006 A/S) @)e eve OOO00S OO0O00E Glalale len! = 


Bry wi ge o 


Sa Oewe 
SuUNSseaug way * 
Mmm Ct'tPt = YORtY 
OOO00V 
OO00S9 
OO008 





Haag! 


Cie 7 Nee 


43 





Pare reducing the flow rate ‘through the tube 
lowers the heat flux and increases the inside convectiv 
thermal resistance. SmicS "enes «GoudensSsing CoetLEicisnt is 
inversely proportional to the heat flux raised to the one- 
third power, the outside condensing resistance décreases. 
Because of this, the inside resistance becomes the dominant 
Bie or in the overall thermal resistance of the tubs. Any 
error in measuring the inside heat transfer ccefficitent will 
greatly amplify the computed error in the condensing reésis- 
tance. Thus, the condenSing coefficient is highly sensitive 
memecne aecciracy of the inside coefficient, particularly at 
low water velocities (and corresponding heat fluxes). 

second, “he test apparatus was subject to radio- 
frequency (RF) interrerence ftcm the power supply. As the 
power was cut back +o lower the heat flux, the silicon 
Penwrolled rectifier chovoped the input signal proportionally 
and emitted the GWCERedePicrE1On 25. RF energy. Despite 
Mmeeores tO shield the cabling to the data-acquisition 
system, it still captured the RF signal with unacceptable 
results, and the thermopile was experiencing instantaneous 
fluctuations of several hundred microvolts. Increasing the 
power supply, however, decreased ths errant signal until, art 
the maximum power input used for th2 high heat flux measure- 
ments, the thermopile showed @&@ near-steady temperatur? 
difference that was within 0.1 °C of the quartz thermometer. 

For these reasons, the Sieder-Tate coefficients 
obtained from the low heat-flux runs were discardéd. 
Averaging the coefficients found from both measurements for 
the high aeat flux case yielded a coefficient of 0.034 + 
0.001. 
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Be. FINNED TUBES 


Figure 5.5 shows a relative DLO. moO rms OUT Side 
Ss 


Sh 
SGendensi:ng coefficient for all six finned <ube compared 
with smooth tube results and the Nusselt lines. Figure 5.6 
shows the same results corrected for area ases due to 
finning. These data were computed usiag the leading cosffi- 
Cient on the inside of 0.034 obtained in the 

section. 


All tubes tested show good enhancement over the smooth 


eaD2. Notice the increase of thea smooth tube over the 
Nusselt line - a result of the inherent vapor shear in the 
mest SeCtion. The increased scatt2r in the lower heat-flux 


range is again a function of the increased Significance of 
errors in measuring the inside thermal resistances. The 
plotting program used neglected ordinate values below zero 
Semeagrea~er chan 105, so the scattering in this pcrtion of 
the plot is actually worse than it appears. Figure 5.7 is a 
foe Of the values oktained for th2 outside heat transfer 
coefficients of the six finned tubes for a constant heat 
flux of 250,000 W/m? (79,000 Btu/hr-tt2) and a pressure of 
88 mm Hg (1.7 pSia). This plot more clearly shows the 
Secemum pitch of 2.5 mn. 

The heat-transfer characteristics for any tube will be 
Benanced by the addition of fins. In the case of purely 
Bonvective heat transfer, the enhancemen* is a function of 
the increased surface area exposed to the fluid nediun. 

Filmwise condensation, however, is dissimilar in «hat 
ete DUIId-up or a ccndensate filn acts as an additional 
thermal resistance for the heat-transfer process. The objec- 
tive of tube enhancement, therefore, is te decrease this 
mem thickness. The surface tension forces in the conden- 
Sate tend to draw the liguid toward the fins, leaving the 
tube surface with a thinner film. The thinner film results 
in a higher heat-transfer coefficient. [Ref. 12]. 
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The optimum pitch found ian this research increases the 
surface of the condenser tube 88%, Ducewecenacising heat- 
transfer coefficients for that pitch were enhanced by Ss 
much as 330%. 

The optimum pitch serves as the right trade-off between 
the attractive forces of the fins on the condensate and the 
channel aréa between the fins to drain the condensates fron 
the tube. A pitch smaller than th? optimum has too narrow a 
gap between fins to efficiently allow condensate run-off. 
This is because the condensate drawn +o adjacent fins only 


BemoiRes to create a thicker film as shown in Figure 5.8. 
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2.0 Sketch of the effect of finning on condensation 
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A. CONCLUSIONS 


1. The test apparatus successfully operates under 
Vacuum conditions with no degrading effects caused by the 
presence of noncondensable gases. 

2- The Sieder-Tate coefficient was deternined +9 be 
0.034 + 0.001 using both an instrumented smooth tube and the 
modified Wilson plot technigue. 

Pe sneeeooc=mun f2nemseen Of 2.5 mm was found for the 
finned tubes tested. 

GW. Heat-transfer enhancements up +0 300% were realized 
while the increase in the tube outer surface area (due to 


mematng) was only 88% over the smooth tube. 


Be RECOMMENDATIONS 


1. Install a second, more reliabl2 pressure transducer 
to replace the manometeér. 

2. Install a larger cooling-watéer sump to eliminate use 
Of the once-through system from the fresh-water tap. 

BepecOlilect Gata USing anh insSr= ticr the test «cube to 
decreass the internal convective rasistance and, thereby, 
decrease the uncertainty 1n the calculated condensing coef- 
mecient. 

4. Take data for the finned tubes while varying fin 
height and thickness (as w21l as pitch). 
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APPENDIX A 


THERMOCOU PLE CALIBRATION 


we GOULPMENT USED 

ae Thermocouple Wire 

Copper-constantan, 0.245-mm (0.01-in) Teflon-coated wire 
was used for all thermocouples. 

be, Calibration bath 

A Rosemont Engineering Model 913A calibration bath was 
used. A schematic representation 9f the bath is shown in 
PeGure A.1. 

Nye Heating:: Blectrical 

we COClang: La quidmi2a trogen 

Notes Once a desired temperature is reached, the temp- 
erature is held constant by rapid cycling between heating 
and cooling. The bath is rated for temperature fluctuations 
@eress than 0.01°C. 

ce. Thermocouple readout 

An HP 3054A Data AcquisitionysControl System was used +o 
Obtain data. Resoluticn of the acquisition systam was 1aV. 

d. Fath temperature measurement 

A platinum rasistance thermometer with an accuracy of 


0.01 CC was used. 


Be PREPARATION 

@. Procedure 

The instrumented tube (with the wall thermccouples 
installed) and the steam thermocouples were immersed in the 
bath as well as the prebes for the guartz thermometer. 

b. Analysis 

The computer program TCAL was used to monitor and store 
all thermocouple readings ona disk. A listing of the 
program is located in Appendix E. 
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36 CALIBRATION PROCEDURE 

a. The bath temperature was sat at about 10°C. 

b. When the bath temperature reached st2ady state, its 
value was entered into the computer. 

ce. The computer automatically recorded and printed all 
thermocouple readings. 

d. The bath temperature was raised in incraments of 10 
C to 90 °C and steps b and c were repeated for each incre- 


Ment e 


4. CALIBRATION CURVES 
ae A least-squares method was used to generate a poly- 


nomial of the form: 
Dr = a, + a,T + a,T2 


wheres Ds one differance between the bath and thermo- 


couple temperatures, and 


1 is the value of the thermocouples obtained using the 
seventh-order polynomial fitted for the Type T thermocourle 
Wire used. Milo meociynOlladl as Listed in the beginning of 
*he program TCAL. 


b. Coefficient values obtained were: 


ag = 4.7338 x 10-3 
a, = 7.6928 x 10-3 
apes soe01s/9 x 10-5 


Sie MepbormoL che Carves 2S £f5und in Figure A.1. The 
curve tor the wali thermocouples giv2s a different reading 
due to the thermal conduction of th2 air temperature through 
the tube. The data used for calibration did not include 
this average value. Only the value obtained for the two 


steam thermocouples was used for all of the thermocouples. 
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aa aS 2 


SAMPLE CALCULATIONS 


fe HEAT LOSS 


A set of calculations was performed to estimate the hear 


joss of the apparatus due +o natural convection. These 
calculations were performed for two main sections: the 
Perler (with and without insulation) Sladen es pa paag tO the 


test section. The minimum power inpat used during the thesis 
was 17 kW. 


ae Boiler (without insulation) 


t4) 
ce 
2 
iD 


1) Givan the following dimensions and properties o 


holler section: 


Biveg = 0.1524 Oo 
Doeg = 9.1567 o 
Lz = 0.4064 n 
Ko = 1.4 W/n-K 


2) Given the following temperatures: 


aoe = 100°C 


Te.g = 95 ¢ 
fee = 25°C 
Meee = (95425) 72 = 60°C 


3) Given the following properties of air at the film <emp- 


ema cure ; 


Ka = 28.74x10-3 W/m-K 


5) 5) 





4) 


Pr = 0.702 

tame= loeZ2ixto-* m-/s 

ome — 05003 K-? 

Compute the Grashof Number for the boiler assuming 


laminar conditions (Ra < 109): 


3) 


6) 


Grg = 9G(T, g ~T.) 13/7? 


3.7h 108 


GIlg 


Compute the Rayleigh Number for th2 boiler: 


Ra, = Gr PI 
Ra - Zeon 10% (< 109) 


Det 2rmine the applicability of a flat-plate analysis 


fescerdaing to che method of Sparrow and Gregg [Ref. 13] 


i7Eo6 < (0.15/18) Grt7* 
2-60 < 7.36 


Compute the average Nusselt Number (via flat-plate anal- 


8) 

ysis): 
Nu = [0.825 + 0.387Ra, t/8/[ 1+(0.492/Pr) 9/16 18/27 J2 
Nu = 81.3 

9) Compute the external heat-transfer coefficient: 
fee = Na k /L 
ho = 5.75 W/m2-K 

10) Compute the external convective thermal resistance: 
R = 1/ (hg 277, 3L3) 
R = 0.436 K/W 

11) Calculate the wall resistance: 
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2) 


Rw 


Rw = 


Neglecting the wall resistance, 


Emel a7 ye) eho lg 


0.007 K/W 


from the boiler: 


Q2g > 


Qs 


(Tie-T 7 (R +R) 


169.3 W 


Bemler "(With douSulat2on } 


95, = 29.24 
Piping 
Q,= 65.0 


Total (worst case) 


Q 


Q 


oe Cor 


235463 0 


(<< 17000 W of 
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APPENDIX C 
STEAM-CLEANING PROCEDURE 


Note: The procedures listed h2re provid? an excellent 
method of cleaning the apparatus, but will cause the walls 
of «he apparatus to heat to temperatures over 100 C (212 F). 
Once this happens, several hours are required before the 
walls can cool down such that they won't superheat the stsam 
generated during an experimental run. Therefore, use this 
method only prior to operating at atmospheric pressure or 
when the system is badly contaminated. 

1. Ensure the water level in the boiler is a+ least six 
inches above the upper ends of the heaters. 

2. Remove the thermopile from the inlet side of the 
test section. 

3. Energize the power supply and adjust the rheostat 
until the voltmeter reads about 100 V. 

4, Close the recirculating valve. 

5. Open the fillydrain valve. 

6. Once vapor begins exiting the fill/drain line, allow 
the system to steam fer several minutes. 

7. j%qTO maximize the steam flow through the tast section, 
fully open the throttle to the auxiliary condenser. 

8. When done, simultaneously close the drain/fill valve 
and open the recirculating valve. 

9. Circulate cooling water through the test tube and 


check for the presence of any dropwise condensation. 


10. Reinstall the thermopile. 
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SYSTEM Seat -U P 


To start the system: 

1. Fill the boiler to at least six inches above the 
upper level of the heaters. 

2. Energize the air ejector for ten minutes. 

3. Circulate water through th2 cooling water sump by 
opening the inlet valve and activating the Siphon. 

4. Energize the power to tha boiler and adjust*+ the 
theostat until the vcltmeter reads 90 ¥V for vacuum opera- 
tion, or 170 V for atmospheric pressures. 

Se ENe=GLze —emecrtoulating pumps and adjust the 
throttl= for the desired flow through the tube. 

6. Open the valve to the auxiliary condenser and adjust 
the flow: about 15% for vacuum operation or 30% for atnos- 
beezic runs. 

7. Observe the steam temperature indicated on th efront 
of the data acquisition unit until the voltage is reached 
corresponding to the saturation temperature of the the 
desired operating pressure. 

8. Energize the air ajector for two minutes. 

9. Adjust the the flow rate through thea auxiliary 
condenser tc bring the apparatus to the desired operating 


pressure. 


To secure the systen: 
1. Secure cooling water to the test tube. 
- Secures power to the boiler. 


iy NBN 


- Open the vent valve. 


B 


-» oecure the water supply to the auxiliary condenser. 
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COMPUTER PROGRAM LISTINGS 


The following pages contain the computer program list- 


ings used tor this thesis: 
SIEDER (page 61) 
WILSCN (page 65) 
DR? (page 71) 


eA L. (page 80} 
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PROGRAM SIEDER 


Se a ee ee ee ee ee | 
(ae et te tae 
COWOoOnNNa S&S WN — 
ooooo°0oat0ooe 


1510 
1520 
io 
1540 
1550 
13560 
1570 
1580! 


“TLE NAME: SIEDER 
DISK NUMBER: 12 


Revise: November 29. 1983 

Come /Cc/ CC7) 

DIM Emf(10).Tuw(5) 

DATA 0.1908609!1 .25727 .94369., -767345.8295. 78025595 .81 
DATA -9247486589.6.97688E11,-2.66192E13.3.94078E14 
READ CC(#) 

Kceu=385 

Ore .0127 

Do=.01905 

Dr=.0158/5 

P= 13335 

Li=.060325 

be, 054925 

PRINTER IS 701 

BEEP 

CLEAR 709 

INPUT “ENTER MONTH, DATE AND TIME (MM:DD:HH:MM:SS)",8S 
OuTeUL 709:-°TD" :8s 
Series: ! 

OUTPUT 709:"TD" 

ENTER 709;:AS 

any USING "“{0X,""Month. date and time: "".14A" SAS 
BEE! 

INPUT “ENTER DISK NUMBER" ,Dn 

PRINT 

PRINT USING “10X,°"’NOTE: Program name: SIEDER’’” 
pan USING “16X.""Disk number = "".DD'’:Dn 

INPUT “ENTER INPUT MODE (1+3054A.2=F ILE)" ,Im 

IF Im=1 THEN 

ee 8 

INPUT “GIVE A NAME FOR THE DATA’ FILE’ .D_fileS 
CREATE BDAT D_filesS.10 

SESE 

BEEP 

ee “GIVE THE NAME OF THE DATA FILE’.D_ Files 

INPUT “ENTER THE NUMBER OF RUNS STORED” .Nrun 
ie "16X,""This analysis was performed for file “",10A";:D_fileS 
Beer 

NeUT "GIVE A NAME FOR PLOT DATA FILE” ,PlotS 

INPUT “ENTER OPTION FOR END-FIN EFFECT (1=Y,0=N)" Ife 
IF Ife=0 THEN PRINT USING “16X.""This analysis neglects end-fin effect” 
IF Ife=1 THEN PRINT USING “16X%,""This analysis includes end-fin effect" 
CREATE BDAT PlotsS,.5S 

ASSIGN @File TO D_filesS 

hee @Filep TO Plot 

Sx=0 

Sy=0 

Sxs=Q 

Sxy=0 

IF Im=1 THEN 

READ DATA THROUGH THE DATA ACQUISITION SYSTEM 
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esl 
1600 
1670 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
70 
.1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880! 
1890! 
1900 
1910 
nec0 
1930 
1940 
1350 
1970 
1980 
1990 
200 
203 
2020 
2030 
2040 
2050 
2051 
2060! 
2070 
2080 
2090 
2,00 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 


IF THE INPUT MODE (Im) = 1 
BEEP 


INPUT “ENTER FLOWMETER READING". Fm 
OUTPUT 709:"AR AF20 AL30 VR1" 
FOR I=0 TO 10 

OUTPUT 709:"AS SA” 

IF I>4 THEN 

Se=0 

FOR K=0 TO 19 

ENTER 709°: 

Se=SetE 

NEXT K 

Emf¢(I)=ABS(Se/20) 

BUSe 


ENTER 709:E 

Emf (I) =ABS(E) 

END IF 

NEXT I 

OUTeUT 713: TtR2ZE" 

WAIT 2 

SMTER 713:1tt 

BoPUT 7133" T2R2E" 

WAIT 2 

ENTER 713;172 

Cel 2io; 1 tR26 

WAIT 2 

ENTER 713:Ti2 

T1=¢€7114T12)*.5 

ELSE 

READ DATA FROM A USER-SPECIFIED FILE IF 
INPUT MODE (Im) = 2 

ENTER @File;sEmf(*#) ,TI.T2.Fm 

END IF 

Tavg=(Ti+T2)#,.5 

Twal1=0 

FOR I=5 TQ 10 ? 
Tu( I-5)=FNTIvsv(Emf ¢I)) 
Twall=Twall+Tw(I-5) 

NEXT I 

Twall=Twall/6 

Cow=F NCpw( Tava) 

Rhow=FNRhow( Tavg) e 
Md=5 .00049E-3+6.9861937E-3*F m 
Md=Nd*(1.0365-1 .96644E-3*71+5.252E-6+TI~2)/.995434 
Mf=Md/Rhow 

Vu=ht/¢(PIl*D01 72/4) 
T2c=T2-¢.0138+.001*Vw"2) 
T2c=T2-.004*Vw 2 

Q=Md*Cow*(T2c-T1) 
Dtw=Q*LOG(Do/Di)d/(2*PI*Kcu*L)*.5 
Twall=Twall-Dtw 
Lmtd=(T2c-T1)/LOG((Twall-Tt)/(Twal!-T2c)) 
Kw=FNKu(Tavg) 

Pi=PL*(Di+tDo) 

P2=P{*<(Di+Dr) 

At=(Do-Di)*PI*Do 

A2=(Dr-Did*PI*Dr 
Hi=Q/(PI*Di~+L«Lmtd) 

IF Ife=0 THEN 

Hic=Hi 

GOTO 2300 
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2200 

270 
2220 
2230 
2240 
2250 
2260 
22/70 
2280 
2290 
2300 
2310 


2320 
2330! 
2340 
2350 
2360 
22/0 
‘2380 
2390 
2400! 
2410! 
2420 
2430 
2440 
Nu 
2450 
2460 
2470 
2480 
2490! 
2900! 
2510 
2520 
2530 
2540 
2550 
25690 
2570 
2580 
2599 
2600 
2610! 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2691 
270 
2/720 
2730 


END IF 

MT=CHi#*P1/(Kcu*Al1)) > .5 
M2=(Hi*P2/(Kcu*A2))°> .5 
Fel=FNTanh(M1l*LIid/CMt#Lt) 
Fe2=FNTanh(M2*L2)/(M2*L2) 
Hic=Q/(PI*Di*(L+Li*Fet+tlLo=Fe2)*Lmtd) 
IF ABS(¢(Hi-Hic)/Hic)>.01 THEN 
Hi=(HictHid*.5 

GOTO 2210 

END IF 

PRINT 


PRINT USING "10X,"°"Position number ; 1 2 3 4 5 
g;,.,"" 


PRINT USING “10X."°"Wall temperature (Deg C) : “",6¢DD.DD,1X)": Tw) 
CALCULATE THE NUSSELT NUMBER 

Nu=Hic*Di/Kw 

Muw=F NMuw< Tavg>? 

Re=Rhow*Vw4D i /Muw 

Cf=(Muw/FNMuw( Twall))*.14 

Pru=FNPru( Tavq) 

X=Re” .8*Prw” .3333*CF 

COMPUTE COEFFICIENTS FOR THE LEAST-SQUARES-FIT 

STRAIGHT LINE 

OUTPUT @Filep;:X.Nu 

PRINT USING "10X,.""Twall in Vay X Nie 
PRINT USING “10X,4¢2D.2D, se ZO). ie MSD, D, a OD ieeetual trl). me. Lmtd. vil. x, 


Sx=Sx+X 

Sy=Sy*Nu 

Sxs=SxstX*X 

Sxy=Sxy+X*Nu 

STORE RAW DATA IN A USER-SPECIFIED FILE IF 
INPUT MODE (Im) = 1 

Ir Im=! THEN OUTPUT ®@File;sEmf(*) ,T?,72.Fm 
BEEP 

J=J+1 

IF Im=1 THEN 
INPUT “ARE YOU TAKING MORE DATA (1=YES.0=NO0)?",Go_on 
Nrun=J 

IF Go_on=! THEN 1570 

ELSE 

IF J<Nrun THEN 1570 

END IF 

Ci=Sxy/Sxs 

Ci=(Nrun*Sxy-Sy*Sx)/(Nrun*Sxs-Sx 2) 

Ac=(Sy-Ci*Sx)/Nrun 

PRINT 

av USING "10X."°°Sieder-Tate Coefficient = °",D.4D":;Ci 
RINT 

PRINT USING “10X.""Least-Squares Line:””” 

PRINT USING "12X,""Slope = eee eaves 2Cl 

PRINT USING “12X,""Intercept =, 1D. SDE, :Ac 

PRINT 

IF Im=1 THEN 

Beer 


ao 


PRIN} USING “10X,°°NDIE: “",7ZZ.°” data runs were stored in file "",8A’:Nru 


The Det iles 


ea03 | 


alse 
lke Bao USING “10X,""NOTE: The above analysis was performed for file "".14A" 
»U_tile 
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2740 
E130 
2760 
2770 
2780 
27390 
2800 
2810 
2820 
2830 


END IF 

Boe mS iiGew tee Plot data are Stored in file “°",.14A7:P lots 
ASSIGN @File TQ * 

ASSIGN @®Filep TO *« 

END 

DEF FNRhow(T) 

Ro=1006.35724-T#*« .774489-T#(2.262459E-2-T#3.03304E-4)) 
RETURN Ro 

FNEND 

DEF FNPru(T) 

Pru=FNCpow< T) *FNMuw¢(T) /FNKu (TD 

RETURN Pru 

FNEND 

DEF FNMuw(T>? 

A=247 .8/¢(T+133.15) 

Muw=2.4E-5#107A 

RETURN Mua 

FNEND 

DEF FNKw(T) 

Kuw=.572183504477+1 .52770121209E-3*T 

RETURN Kw 

FNEND 

DEF FNIvsv(Emf) 

cum /Ce/ CC7) 

Sum=C (0) 

FOR I=! TQ 7 

Sum=SumtC (I) «Emf *I 

NEXT [ 

RETURN Sum 

FNEND 

DEF FNCpw(T) . 
Cew=(4.21120858-T#(2.26826E-3-T*¢4.42361E-S+T*2.71428E-7)))=1000 
RETURN Cpw 

FNEND 

DEF FNTanh(X) 

P=EXP(X) , 

G=EXP(-X) 

Tanh=(P+Q)/(P-Q) 

RETURN Tanh 

FNEND 
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PROGRAM WILSON 


1000! FILE NAME: WILSGN 

1010! REVISED: December 5. 1982 

fo20! 

moO COM /Cc/ CC?) 

1040 DATA 0.10086091 .25727 .94369.-767345 .3295, 78025595 .31 
1050 DATA -9247486589.6.97688E11.-2.66192£13.3.94078F14 
1060 READ CC) 

1070 DIM Emf (4) 

mwe0 L=.130175 

1090 Lt=.060325 

1100 L2=.034925 

my10 Do=.01905 

1120 Di=.0127 

mi3o0 «Dr=.015785 

1740 Keu=3385 

1150 Rm=Do*LOG(Do/Di)/(2*Kcu) 

1160 PRINTER IS 701 

1770 BEEP 

‘1180 CLEAR 709 

1190 INPUT "ENTER MONTH. DATE, AND TIME (MM:9D:HH:MM:SS",BS 
1200 OUTPUT 709:"TD":Bs 

1210 Jp=0 

1220 QUTPUT 709;"TD” 

1230 ENTER 709:AS 


1240 PRINT USING “i0X,""Month, date and time : “",14A";:AS 
#250 BEEP 

1260 INPUT "ENTER DISK NUMBER",Dn 

1270 PRINT 

1280 PRINT USING “10X.""NOTE: Program name : M_WILSON"'” 
1290 PRINT USING "16X,""Disk number Se BIR hg 

1300 BEEP 


Pores INPUT “ENTER INPUT MODE (7 =3054A,2=F ILE)” .Im 
1320 IF Im=1 THEN 

fee0 BEEP , 

1340 INPUT “GIVE A NAME FOR THE DATA FILE”’,.D_ files 
1350 Bice |e BDAT D_filesS,i0 


ia70 BEEP 
bas0 INPUT “GIVE THE NAME OF THE DATA FILE”.D files 
1390 eet USING "16X,""This analysis is for data in file "",14A";D_ files 


EEP 
14470 INPUT "ENTER THE NUMBER OF RUNS STORED” .Nrun 
1420 END IF 


1440 INPUT “GIVE A NAME FOR PLUT-DATA FILE”,PlotS 
1460 INPUT "ENTER OPTION ¢(1=QCT,2=T-PILE.3=AVE)”",Itm 


1480 INPUT “ENTER OPTION FGR END-FIN EFFECT (1=Y,0=N)".I fe 

1490 IF Itm=1 THEN PRINT USING "16X,""This analysis uses QCT readings””” 

1500 IF Itm=2 THEN PRINT USING “16X,""This analysis uses T-PILE readings”’” 
1510 IF Itm=3 THEN PRINT USING "16X,°"This analysis uses average of QCT and [-P 
ILE readings””” 

1520 IF Ife=f THEN PRINT USING “16X,""This analysis tncludes end-fin effect"’” 
1530 IF Ife=0 THEN PRINT USING “16X,""This analysis neglects end-fin effect’”’” 
1540 CREATE BDAT PlotS,10 

1550 ASSIGN &Filep TO Plot$S 

1560! Ciu=.040 

1570! Cil=.028 


65 





1580 
1590 
1600 
1510! 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740! 
1750! 
1760 
1770 
1780 
1790 
1800 


Sxy=0 

PRINT 

PRINT USING "10X."“Iteration number 

IF Jy=0 OR Jp=1 THEN 

PRINT 

erin! USING “12xX.°"°T1 T2 iSat 
END IF 

ASSIGN @File TO D_files 

IF Im=1 AND Jy=0 THEN 


= acuity, ©) Bs a ae oe 


imtd Vu XK de tes 


READ DATA THROUGH THE DATA ACQUISITION SYSTEM 


eeu INPUT MODE (Im) = 1 
INPUT “ENTER FLOWMETER READING" ,Fm 
BDUTPUT 709:"AR AF6BO ALB3" 
OUTPUT 709:"AS SA” 

Etp=0 

FOR I=! TO 20 

ENTER 709;3E¢ 

Etp=Etp+Et 

NEXT I 

Etp=Etp/20 

OUTPUT 709:""AS SA” 
Ptran=0 

FOR I=1 TO 50 

ENTERS /09:Pt 
Ptran=PtrantPt 

NEXT I 

Peran=rtran/50 

OUTPUT 709:"AS SA" - 
ENTER 709:Bvol 

GUTPUT 709:"48S SA" 

ENTER 709:Bamp 

OUTPUT 709;"AR AF2ZO AL24" 
FOR I=0 TO 4 

OUTPUT 709;"AS SA” 

ENTER 709:Emf (I> 

Emf (I) =ABS(Emf (I) > 

NEXT I 

OUTPUT 713;3"TIR2ZE” 

WAIT 2 

ENTER 7133711 

DUTPUT 713:"T2R2E” 

WAIT 2 

ENTER 713;T2 

OUTPUT 713;"TIR2ZE” 

WAIT 2 

ENTER 713:T12 
Tiz(T11+T12)*.5 

CLEAR 713 

ELSE 


READ DATA FROM A USER-SPECIFIED FILE IF INPUT MODE (Im) = 2 
ENTER @File;Bvol,Bamp.Ptran,Etp.Emf(=),.Fm,71,12 


END IF 
Tsat=FNIvsv¢ (Emf (0) +Emt (1))>*.5) 
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2780 
27930 


Ti=FNTvusvtEmf(2)) 
Grad=FNGrad((T{+T2)*.5) 
To=Tit+tABS(Etp)/(10*Grad)*1.E+6 
IF Jy=0 THEN 
Erl=ABS(Ti-T1) 
PRINTER IS 1 
Gren) USING ~°"'T 1 
PRINT USING °° °T:i 
IF Erl>.S THEN 
BEEP 

PRINT “QCT AND TC DIFFER MORE THAN 9.5 C” 


at Dees) Dass ia 
ene DO. . Tt 


a 


BEEP 

INPUT “OK TG GO AHEAD (i=Y¥,0=N)?" Ok! 

END IF 

feat! USING “""DT <@CT) Geng co | | 
ene lotnee Di <T-PILE) = "".Z.3D":To-Ti 
tae l=UeAND Eri>.> THEN 3600 

mreznbo  le-11)-CTo-Tid)/(12-T)) 

eeere>. 0a THEN 

BEEP 

PRINT “@CT AND T-PILE DIFFER MORE THAN 52” 


Beer 

INPUT "OK TQ GQ AHEAD (1=¥,0=N)?" ,Ok2 
IF Ok2=0 AND Er2>.05 THEN 3600 

END it 

PRINTER IS 701 

END IF 


PCALCULATE THE LOG-MEAN-TEMPERATURE DIFFERENCE 


IF Itm=1 THEN 

TF=T1 

Tl=T2 

END IF 

IF Itm=2 THEN 

Tf=Ti 

Tl=To 

END IF . 
IF Itm=3 THEN 

TE=HCTI471)*.5 

Tl=(T2+To)d#.5 

END IF 

Tava=(Tf+T1l)*.5 

Trise=l1-Tf 
Lmtd=Trise/LOG((Tsat-Tf)/¢Tsat-Tl))? 
Cow=FNCpwlTavg) 

Rhow=F NRhow( Tavg) 

Kuw=F NKw( Tavg) 

Muwa=F NMuw( Tavg) 

Pru=FNPru< Tavg) 
Mdt=5.00049SE-3+6.9861937E-3*Fm 
Md=Mdt#(1.0365-Tf#¢(1.96644E-3-Tf#5.252E-6))/.995434 
VF=Md/Rhow 

Vw=VF/(PI*Di*2/4) 
Trise=Trise-.004*Vw~2 
Q=Md*Cpw*Trise 

Qp=G0/(PI*Do*L) 

Uo=Op/Lmtd 

Re=Rhow*Vu2D1i/Muva 

Fel=0 

Fe2=0 

Cfr=1 

Two=Tsat-5 
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210 


Tfilm=Tsat/2+Two/2 

KF=FNKu(T film) 

Rhof =FNRhow( Tf ilm) 

Muf =FNMuw(Tfilm) 

Hf qp=FNHfFg(Tsat)+.68"FNCpu( Tfilm)=(Tsat- Two) 
New=Kf*(Rhof ~2*9.799*Hf gp/ (Muff *Do*Qp))* .3333 
Ho=.655*New 

Twoc=Tsat-Qp/Ho 

IF ABS((Twoc-Two)/Twoc)>.001t THEN 
Two=Twoc 

GOTO 2800 

END IF 

Cf=1.0 

Omega=Re~ .S*Pru > .3333«Cf 
Hi=Ku/Di*Ci*Umega 

IF Ife=0 THEN 3040 

P1=PI=*(Di+Do) 

P2=PI*(0.i+Dr) 
At=(Do-Di)*PI*«(Di+Do)*«.5 
A2=(Dr-Did«PI*<(Di+Dr>*.5 
M1=CHi*P1/(Kcu*A1))7.5 
M2=(Hi*P2/(Kcu*A2))°.5 
Fel=FNTanh(Ml*Li)>/CM1#*L1) 
Fe2=FNTanh(M2*L2)/CM2*L2) 
Dt=Q/(PI*Dix*(L+L1*Fel+L2*Fe2) *Hi> 
Cfc=(Muwa/FNMuw(TavgtDt))~.14 

IF ABS((CFKe-CF)/CFC)>.01 THEN 
Cre(CftCfc)*.5 

GOTO 2930 

END IF 

X=Do*New/ (Omega*Kw) 
Y=New*(1/Uo-Rm) 


Sue eecunrUte COEFFICIENTS FOR THE LEAST-SQUARES-FIT STRAIGHT LINE 


3130 
3140 
3150 
3160 
31/0 


IF Jp=1 THEN QUTPUT #@Filep:X.Y 
Sx=Sx+X 

Sy=SytY 

Sxs=SxstX#X 

Sxy=Sxy+X*VY 


31380! STORE RAW DATA IN A USER-SPECIFIED FILE IF INPUT MODE (Im = 1 


3190 IF Im=1 AND Jys=0 THEN OUTPUT ®File:Bvol.Bampo.Ptran.Etp.Emf¢(=) ,Fm.7T1,172 
3200 IF Jy=0 OR Jp=1 THEN PRINT USING °8X.5¢(2X,.3D.DD) .2¢2X,D.SD)"“3Tf,Tl.Tsat.im 
td.,Vw,X.Y 

3210 BEEP 

6220) J=J+i 

3230 IF Im=t AND Jy=0 THEN 

3240 INPUT "DO YOU HAVE MORE DATA (1=Y,0=N)?",Go_on 

3250 Nrun=J 

3260 IF Go_on=1 THEN 1730 

e270 ELSE 

3280 IF J<Nrun THEN 1730 

S290 END IF 

3300 Sl=(Nrun*Sxy-Sy*Sx)/(Nrun*Sxs~Sx °2) 

3310 Ac=(Sy-S1*Sx)/Nriun 

Soa2ueCcic=1/S1i 

Socun J;=Jj+1 

3340 IF Jp=! THEN Jp=2 

3350 IF ABS<(Cic-Ci)/Cic)>.001 THEN 

3360 Ci=(Cict+Ci)*.5 

3370 PRINT USING “10X,.""Intermediate Sieder-Tate coefft = "",Z.4D":Ci 
3380 GOTO 1600 
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SoU EL Sc 

3400 IF Jp=0 THEN Jp=] 
3410 END IF 

3420 IF Jp=! THEN {600 
eas0) Ci1=¢(Ci+Cic)*.5 


3440) PRINT 

3450 PRINT USING “10X,.""Sieder-Tate coefficient ee el 
3460 PRINT 

3470 PRINT USING "“10X."’Least-Squares Line: 

3480 PRINT USING “10X,"" Slope ee) 2 ODE. 3S] 

3490 PRINT USING “10X,"" Intercept = "",MZ.5DE,°:Ac 

3500 PRINT 

3510 IF Im=! THEN 


3o20 BEEP 

eeaue PRINT USING “10X,°°NOTE: “",7ZZ.°° data runs are stored in file “°”,8A°:J.D_ 
files 

3540 ELSE 

3550 PRINT USING "“70X."°°NOTE: Above analysis was performed for data in file "", 
10A";D_files 

*3560 END IF 

3570 PRINT USING “16X.°°Plot data are stored in file “",10A”";:PlotS 
32580 ASSIGN @File TO * 

35390 ASSIGN &Fileno TO « 

3600 END 

3610 DEF FNRhow(T) 

3620 Ro=1006.35724-T#¢.774489-1T*¢(2.262459E-2-1T+#3.03304E-4)) 
3630 RETURN Ro 

3640 FNEND 

3650 DEF FNPru(T) 

3660 Pruw=FNCpow(T)*FNMuw?( T)/FNKuCT) 

3670 RETURN Pru 

3680 FNEND 

3690 DEF FNMuww(T) 

3700 A=247 .8/¢(T+133.15) 

3710 Mu=2.4E-5*10°A 

3720 RETURN Mu = 

3730 FNEND 

3740 DEF FNKuw(T) 

Byoo X=(7+273.15)/273.15 

3760 Kw=-.92247+X*(2.8395-X*¢(1.8007-X*¢ .52577-.07344*X))» 
3770 RETURN Kw 

3780 FNEND 

3790 DEF FNIvsv(Emf) 

Bouu COM /Ce/ Cc?) 

3810 Sum=C(0) 

See0 FOR [=! 10 7 

3830 Sum=Sum+C¢(I)«Emf*I 

3840 NEXT [ 

3850 RETURN Sum 

3860 FNEND 

3870 DEF FNCpw(T) 

3880 Cpw=(4.21120858-T*(2,.26826E-3-T*(4.42361E-5+2.71428E-7)))*1000 
3890 RETURN Cpu 

3900 FNEND 

3910 DEF FNTanh(X? 

3920 P=EXP(X) 

3930 GQ=EXP(-X) 

3940 Tanh=(P+Q)/(P-Q) 

3950 RETURN Tanh 

3960 FNEND 

3370 DEF FNGrad(T) 
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3980 
5298 
4000 
4010 
4020 
4030 
4040 
4050 
4050 
4070 
4080 
4990 


CU Cec, ) 

Grad=37 .9853+.104388*T 
RETURN Grad 

FNEND 

DEF FNFvst(T) 
F=466.444+7T*#¢(7.09451-T#1t .65808E-2) 
RETURN F 

FNEND 

DEF FNHfg(T) 

Hf g=2477200-2450*(T-10) 
RETURN Hfg 

FNEND 
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PROGRAM DRP 


Jh—= OW ONO ON 
Ceoocn”rooodoqooe& 


C 


t FILE NAME: DRP 
' REVISED: November 18. 1983 


Gam /Ce/ CC7) 
Die Em? ©7109 
Vaan ONG50G!l,.25/27 .94369.-767345.38295. 738025595.81 


DATA -9247486589,6.97688E+11,-2.66192E+13,3.94078E+14 

READ C(*) 

Di=.0127 t Inside diameter of test tube 

Do=.01905 ' Qutside diameter of test tube 

Ure=.015375 ' Qutside diameter of the outlet end 

Dssp=.1524 ' Inside diameter of stainless steel test section 
Ax=PI*Dssp “2/4-P!1«Do#L 

P=. 1301/5 ' Condensing length 

Li=.060325 ' Inlet end "fin length” 

2=.034925 ‘woutiet ende. tin length 

Kcu=385 ' Thermal conductivity of Copper 

Ci=.034 ' Sieder-Tate coefficient 
Rm=Do*LOG(Do/Di)d/(2*Kcu) '! Wall resistance based on outside area 
PRINTER IS 701 

CLEAR 709 

BEEP 

INPUT “ENTER MONTH, DATE AND TIME (MM:DD:HH:MM:SS)”,Dates 
GUTPUT 709:"TD”’;:Dates 


Heme A093 TD 
ENTER 709;Dates 


PRINT Hontnmeagate ana time : ‘Dated 
PRINT 

PRINT USING "10X.”°’°NOTE: Program name : DRP*’” 
BEEP 

INPUT "ENTER DISK NUMBER’,Dn é 

BaTNT USING “16X,°°Disk number = *",DD’:Dn 

BE 

INPUT “ENTER INPUT MODE (1=3054A,2=F ILE)’, Im 
IF Im=1 THEN 

BEEP 

INPUT “GIVE A NAME FOR THE RAW DATA FILE” ,D_ files 


CREATE BDAT DufileS.15S 
ASSIGN aFile TO D_files 
Bee P 


INPUT “ENTER GEOMETRY CODE Ci=r INNED.O=PLAIN)"”,Ifq 
OUTPUT @File:Ifg 
IF Ifg=0 THEN 


BEEP 
INPUT “WALL TEMPERATURE MEASUREMENT (1=Y,0=N)", cut 
BSE 
Been 
2 ae Pee iChowtOiH AWD HEIGHT ,Fp,Fu.rh 


IF Ifg=0 THEN OUTPUT $File;Iut 


IF Ifg=1 THEN OUTPUT @File:Fp,Fu.Fh 

BESE 

BEEP 

INPUT “GIVE THE NAME OF THE EXISTING DATA FILE’ ,D_fF:les 

PRINT USING "16X,"""This analysis was performed for data in file "".,.10A":D 


1) 





TlLLied 

1550 “BEEP 

1S60 INPUT “ENTER THE NUMBER OF RUNS STORED". Nrun 

1570 ASSIGN @File TO D_files 

1580 ENTER @File:Ifg 

1990 IF Ifg=0 THEN ENTER @File:Iut 

1600 IF Ifg=! THEN ENTER aFile:Fp,Fu.Fh 

1610 END IF 

1620 IF Ifg=0 THEN 

1630 BEEP 

1640 INPUT “WANT TO CREATE A FILE FOR Nr vs F (1=#Y,02N)9" Inf 

ti50 ELSE 

1660 Inf=0 

1670 END IF 

1680 IF Inf=! THEN 

1690 BEEP 

1700 INPUT “GIVE A NAME FOR Nr vs F FILE” Nri#S 

mau CREATE BDAT Nr#$S.2 

1720 ASSIGN O@Nrf TO NrfS 

m7 30 8 68END IF 

1740 BEEP 

1750 INPUT "ENTER OPTION (1=QCT,2=T-PILE.3=AVE)",Itm 

1760 BEEP 

1770 INPUT “ENTER OPTION FOR END-FIN EFFECT ¢(!=Y¥,.0=N)” Ife 

1780 IF Itm=1 THEN PRINT USING "16X%,""This analysis uses QCT readings”’” 

1790 IF Itm=2 THEN PRINT USING “!6X,""This analysis uses T-PILE readings”’”" 
1800 IF Itm=3 THEN PRINT USING "16X.""This analysis uses average of QCT and IT-P 
ILE readings"’*” 

1810 IF Ife=! THEN PRINT USING “16X%.°"°This analysis includes end-fin effect” 
1820 IF Ife2M THEN PRINT USING "16X.°"%This analysis gee nes end-fin effect’’" 


1830 PRINT USING "16X,""Sieder-Tate coefficient = ars Or 
1840 BEEP 

1850 INPUT “GIVE A NAME FOR PLOT DATA FILE”,P_fileS 
1860 CREATE BDAT P_files,5 

1870 ASSIGN 8Filep TO P_filesS , 

1880 IF Iwt=1 THEN 

1890 BEEP 

1900 INPUT "GIVE A NAME FOR WALL TEMPERATURE FILE’ Wtfs 
1910 CREATE BDAT Wtfs,.5 

1920 ASSIGN Filet TO Wtfs 

1930 END IF 

ou) BEEP 

1950 INPUT “ENTER OUTPUT VERSION (1=SHORT,2=LONG)”. Tov 
1960! IF Im=1 THEN 

fo 0! QUTIPUT *®File:Ifa 

1980! IF Ifg=0 THEN QUTPUT ®@File;Iuwt 

1990! IF I[fgq=! THEN CUTPUT 3File:Fp.Fu.Fh 

2000! ELSE 

eo! ENTER ®File:Ifq 

2020! IF Ifg=0 THEN ENTER ®@File;Iuwt 

POgOt iF Iifg=! THEN ENTER ®File:Fp,Fu.Fh 

2040! END IF 

2050 IF Ifg=9 THEN 

2060 PRINT USING “16X,"°"Tube type peReAIh = 

2070 ELSE 

2080 PRINT USING “16K,""“Tube type ; RINNED 
SeeveeeRIN?T USING “16X.”°"Fin pitch, width. and height (mm): “",DD.D,2X.Z.DD.2xX,.Z 
Pe cFO.fw.F A 

2100 END IF 

evi )«=6J=0 


ete) IF lov=} THEN 


We 





2130 
2140 
=a 
E50 
2170 
2130 
2130 


2200 
2210 
2220 
22230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2520 
2330 
2340 
2350 
2360 
2370 
23580 
2390 


PRINT 
Pre iat—) o1REN 
PRINT USEING “10X.°" Data Vu Uo Ho Qp 


PRINT USING “10K,"°" # (m/s) (W/m 2-H) (W/m 2-K) (W/m *2) 
BaSe 


PRINT USING "10X.""Data Vu Uo Ho 
PRINT USING "10X,"°" # (m/s) (W/m 2-K) (W/m 2-K) 


END IF 
END IF 
Sx=0 
Sy=0 
Sxs=0 
Sxy=90 


Repeat:! 


J=J+1 

IF Im=! THEN 

BEEP 

INPUT “LIKE TO CHECK NG CONCENTRATION CI=Y,0=N)?" Na 
BEE 


P 
INPUT “ENTER FLOWMETER READING’ .Fm 
OUTPUT 709:"AR AF&O AL63 VRS" 
OUTPUT 7039:"AS SA” 
ENTER 709:Ete 
GUTPUT 709:"AS SA” 
Vtran=0 
FOR I=! TQ 50 
ENTER 709;:Vt 
Vtran=VtrantVt 
NEXT I 
Vtran=Vtran/50 
Guleut 705: AaS SA” 
ENTER 709;Bvol 
OUTPUT 709:"AS SA” 
ENTER 709;Bamp ao 
IF Iwt=0 THEN OUTPUT 709:"AR AF2Z0 AI24 VRI" 
IF Iwt=1 THEN QUTPUT 709;"AR AF20 AL30 VRI" 
IF Iwt=0 THEN Nn=4 
IF Iwt=1 THEN Nn=10 
FOR I=0 TO Nn 
OUTPUT 709:"AS SA” 
IF I>4 THEN 
Se=0 
FOR K=0 TO 10 
ENTER 709;E 
Se=SetE 
NEXT K 
pee oer lO) 


ESE 
BNIER 7O95E 
Emf (I) =ABS(E) 
END IF 


NEXT I 

OUTPUT 713; "TIRZE” 
WAIT 2 

ENTER 713:711 
OUTPUT 713:"T2R2E" 
WAIT 2 

ENTER 713;72 
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Vy 


(meso 


Qe 
(W/m 2) 


Nr 


Wig 
Cn sd 





Panu oWUteol 7133 °TIRZE” 
2720 WAIT 2 
2730 ENTER 713:TI2 
eyeg) T1=C(711+T129*.5 
2750 IF Ng=0 THEN 2800 
2760 BEEP 
2770 INPUT “ENTER MANOMETER READING (CHL,HR.HRWD” HL Hr.Hry 
2/780 Phg=H!i+Hr 
2790 Pwater=Hr-Hru 
2800 ELSE 
2810 IF Ifg=! OR Iwt=0 THEN 
2820 ENTER @File:Bvol,Bamp.Vtran.Etp.Emf (0) .Emf Cl) .Emf (2) ,Emf (63) .Emf (4) Fm. T1.T 
2.Phg,Pwater 
2830 END IF 
2840 IF Ifg=0 AND Iwt=1 THEN ENTER @®File:Bvol.Bamp.Vtran.Etp.Emf(*) Fm. Tl.t2.Ph 
g,Pwater 
2850 IF J=! GR J=!0 OR J=20 OR J=Nrun THEN 
2860 Ng=1 
mero 60 LSE 
2880 Ng=0 
2890 END IF 
2900 END IF 
2910 Tsteam=FNIvsv((Emf(0)+Emf(1))«.5) ' COMPUTE STEAM TEMPERATURE 
2920) Troom=FNIvsviEmf (3)) 
2930 IF Iwt=! THEN 
2940 Twm=0. 
Zoo0 FOR [#0 TO 5S 
2960) Tw l)=FNIveviEmf(I+5)) 
2970) Twm=TumtTu( 1) 
2550 NEXT I 
2990) 3 Tum=Tum/6 
3000 END IF 
3010 Tcon=FNTvsviEmf (4)) 
3020 Psat=FNPvst(Tsteam) 
3030 Rohg=13529-122*(Troom-26.85)/50 
3040 Rowater=FNRhow<( Troom) 
3050 Ptest=(Phg*Rohg-Pwater*Rowater)*9.799/1000 
3060 Pmm=Ptest/133.322 
3070 Pkm=Ptest*1.£-3 
S060 Pks=Psat*1.£-3 
3090 Pkt=FNPysv(Vtran)*#!.£-3 
3100 Tsat=FNIvsp(Ptest) 
3110 Vst=FNUvst(Tsteam) 
ome 0) Peng=(Ptest-Psat)/Ptest 
onoU) Ppst={-Ppng 
3140 Mfng=1/¢1+18.015/28.97<Psat/(Ptest-Psat)) 
3150 Vfng=Mfng/(1.608- .608*Mfnq) 
3160 Mfnag=Mfng*!100 
3170 Vfng=Vfng*100 
BEEP 


3190 IF Tov=2 THEN 

3200 PRINT 

3210 PRINT USING “10X,"““Data set number = OD es 

3220 PRINT 

3230 END IF 

3240 IF Tov=2 AND Naq=I THEN 

pay RINT WSTNG  “10x,°" 2 Paat Ptran Tmeas Tsat N 

3260 PRINT USING "10X."" (mm) (cea (kPa) (k Pad (C) (C) Molal 
Ss 


Severe Rii USING “i0X,5(3D.9D.2X) .2¢3D.0D.2X).2(N3D.D.2X)" :Pmm.Pkm.Pks.Pkt.Iste 
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an, Tsat.Vfng.Mfne 

3280 PRINT 

S290 END IF 

So00 IF Mfng>.5 THEN 

3310 BEEP 

3320 PRINT 

3330 IF Im=1 THEN 

3340 BEEP 

3350 PRINT 

3360 PRINT USING “t0X,.""Energize the vacuum system """ 
3370 BEEP 

3380 INPUT “GK TO ACCEPT THIS RUN (13=Y¥,0=N)?",Ok 
3390 IF Ok=0 THEN 

3490 BEEP 

3410 DISP “NOTE: THIS DATA SET WILL BE DISCARDED!! * 
3420 WAIT 5S 

3430 GOTO 2280 

3440 END IF 

3450 END IF 

3460 END IF 

3470 IF Im=! THEN 

3480 IF Ifgq=1 OR Iwt=0 THEN 

3490 OUTPUT &File:Bvol,Banmp,.Vtran.Etp.Emf (0) .Emf (1) Emf <2) .Emf (3) ,Emf (4) .Fo,TI, 
T2.Phg.Puwater 

3500 END IF 

3510 IF I#g=0 AND Iuwt=! THEN OUTPUT 3Firle:Bvo! ,Bamp.Vtran.Etp.Emf(*) .Fm.T1.72.P 
hq ,Puwater 

3520 END IF 

3530 IF Ifq=0 AND Iwt=! THEN OUTPUT aFilel:Tu(*) 
3540! ANALYSIS BEGINS 

3550 Ti=FNTvsv¢Emf (2) ) 

3560 Grad=FNGrad((T1+T2)*.5) 

3570 =To=Tit+ABS(Etp)/(10*Grad) #1 .E+6 

3580 Erl=ABS(Ti-T1) 

35990 PRINTER IS 1 


Sau PRINT USING “""T1 «(GCT) — D0eeD it 
eeu PRINT USING """T1 (TCO) ee Oa Vt 
See IF Eri>.S THEN 

Seo BEEP 

Seaver niINT “GCT AND TC DIFFER BY MORE THAN 0.5 C™ 
Bo50) BEEP 

3660 INPUT “OQK TO GO AHEAD (1=Y.0=N)?" 0k! 

3670 END IF 

3680 PRINT USING "““"“DT (GCT) meee) boo ih 


Saar enin! USING “""DT (T-PILE) = “",Z.3D"%:To-Ti 

3700 IF Ok!=0 AND Ert>.S THEN 5100 

emo, Er2-ABS(<(12-11)-¢To-Ti))/¢(T2-T1) 

Syece IF Er2>.05 THEN 

you §«6BEEP 

3740 aw SQGimeadp Ver iLe DIFFER BY MORE THAN 52%" 
Pp 


3760 INPUT “OK TO GO AHEAD ¢1=Y,0=N)?" ,Ok2 
3770 IF Ok2=0 AND Er2>.05 THEN 5100 

3780 END IF 

3790 PRINTER IS 701 

3800 IF Itm=t THEN 

3810 Tli=T1 

3820 T20=T2 

3830 END IF 

3840 IF Itm=2 THEN 

3850 Tli=Ti 
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4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
4460 


T2o02Ta 

END if 

te eitm=3 THEN 

P= 1+ lr). 5 

T20=(T2+To)*.5 

END IF 

Tavea=<(liitlzo)*.5 

Cow=—FNCow( Tava) 

Rhow=F NRhow( Tavq) 

Md=5 .00049E-3+6 .9861937E-3*Fim 
Md=Md*(1.0365-1.96644E-3#Tavgt5.252E-6*Tavg> 2)/.995434 
Mf=Md/Rhow 

Vw=hf/(PI*Di~2/4) 

T2o0= F20-¢€.0138+.001*V¥w 2) 
Q=Md*Cpw*(TZ0-Tti? 

Qp=Q/(PI*Do*L) 

Kw=FNKuw(Tavg) 

Muw=F NMuw( Tavg) 

Re1=Rhow*Vw*D i / Muu 

Prw=FNPru(Tavg) 

Fel=0. 

Fe2=0. 

Cfr=1. 
Hi=Ku*Ci/Di*xRer* .S*Pru .3333*CF 
Dt=Q/(PI*Dix(L+L1*Fel+L2*Fe2)*Hi) 
Cfco=(Muw/FNMuw( Tavg+Dt))>.14 

IF ABS((CFC-CFIO/CFc)>.01 THEN 
Ch=(CF+C fod *.5 

GOTO 4090 

END Te 

IF Ife=) THEN GOTO 4250 
P1=PI*(D1+Do) 
A1=(Do-Di)*PI*(Di+Do)*.5 
Mi=CHi*P1/(Kcu*A1))*.5 
P2=PI*(Di+Dr)d 
A2=(Dr-Di)*PI*(DitDr)*.5 y 
M2=(Hi*P2/ (Kcu*A2))°.5 
Fel=FNTanh(M1*L1)/(41*L1) 
Fe2=FNTann(M2*L2)/(M2*L2) 
Lmtd=(T20-T1i)/LOG((Tsteam-T11)/(Tsteam-T2o) ) 
Yo=O/(Lmatd*PI*Do#*L) 
Ho=!/¢C1t/Uo-Do#*L/(Di*(Lt+Li *Fel+L2*Fe2) *Hid-Rm) 
Dtc=Q/(PI*Dix¢(L+L1*Fel+lL2*Fe2) *H1) 
TF ABSC(Dtc-Dt)/Dtc)>.01 THEN 4090 
Hf g=FNHfg(Tsteam) 

Two=Tsteam-Qp/Ho 
Tfilm=Tsteam/3+2*Two/3 

KF=FNKu(T Film) 

Rhof =FNRhow(TFilm) 

Muf =FNMuw( Tfilm) : 
Hpq=.651*KfF*(Rhof “2*9.81*Hfg/ (Muff *Do#*Qp)) ~.3333 
Y=Hpq*Qp* .3333 

X=Qp 

Sx=Sx+X 

Sy=Sy+Y 

Sxs=Sxst+X*2 

Sxy=Sxytk*Y 

OUTPUT &Filep:Qp.Ho 

Q@li=500 ! TO BE MODIFIED 
Qloss=Q1/(100-25)*#(Tsteam-Troom) ! TG BE MODIFIED 
Hf c= NHF (Tcon) 
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4470 Hf=FNHF(Tsteam) 

4480 Mdv=0 

4490! Bp=(Bvol<100)°2/5.75 ! BOILER POWER I[N Watts 
4500 Bp=(Bvol#100)°2/5.76 

4510 Mdvc=((Bp-Qloss)-Mdv*(Hf-Hfe))/Hfa 

4520 IF ABS¢(Mdv-Mdvc)/Mdvc)>.01 THEN 

4530 Mdv=(Mdvt+Mdvc)*.5 

4540 GOTO 4510 

4550 END IF 

4560 Mdv=(Mdvt+Mdve)*.5 

4570 Va=FNVvst(Tsteam) 

4580 Vv=Mdv#Vg/Ax 

4590 IF Inf=! THEN 

4600 F=(9.799*Do«*Muf«Hfg)/ (Vv “2*KF*(Tsteam- Two) ) 
4610 Nu=Ho«Do/Kf 

4620 Ret=Vv*Rhof =Do/Muf 

4630 Nr=Nu/Ret*.5 

4640 END IF 

4650 IF Inf=! THEN QUTPUT @Nrf:F.Nr 

4660 IF lIov=2 THEN 


4670 PRINT USING “10X."" T (Inlet) Delta=1 

4680 PRINT USING “10X,"" QCT iG OCr eer tba 

eee in uolNG “10X,.2¢DD.0D,2X%) .2¢2.30.2X)" 3 01,7T1.7T2-T1.fo-T1 

4700 ce USING “10X,"" Vu Rei Hi Uo Ho q 
V 


See PRINT USING ~10X.Z.D0D.1X.5<¢NZ.3DE.1X).MZ.D0D° :Vu.Rei .Hi,Uo.Ho.Qp,Vv 
4720 END IF 

4730 IF Iov=! THEN 

4740 IF Inf=i THEN 

Pee ee MW oING T1X,0D.2%,2.DD.2X,2(90.D.2X).Z.3DE.1X.Z.DD,.2¢1X.3D.D0D)° 35. Vu. 
Uo.Ho.Qp.Vv.F .Nr 

arou, ELSE 

Seder aiNt USING “11X,DD.2X.Z.DD.2X.2(ND.4DE.2X) .Z2.3DE.3X,Z.0D":J,Vu,Uo.Ho.9p.V 


V 

4780 END IF 

4790 END IF , 

4800 IF Im=! THEN 

43810 BEEP 

4820 INPUT “WILL THERE BE ANUTHER RUN (1=Y.0=N)?°,Go_on 

4830 Nrun=J 

4840 IF Go_on=! THEN Repeat 

485 ELSE 

4860 IF J<Nrun THEN Repeat 

4870 END IF 

4880 IF Ifg=0 THEN 

4890 PRINT 

4900 Sl=(Nrun#*Sxy-Sy*Sx)/(Nrun*Sxs-Sx ~2) 

4910 Ac=(Sy-S1*Sx)/Nrun 

4920 PRINT USING “10X.”""Least-Squares Line for Hnu vs gq curve:'”” 

4930 PRINT USING "10X,""" Slope 20 MD. 40E +S] 

4940 PRINT USING “"10X,"" Intercept = "",MD.4DE":Ac 

4950 END IF 

43960 IF Im=! THEN 

4970 BEEP 

4980 PRINT 

eae aN Bolo 104, «nOre. 2.22,  °£4xdata runs were stored in file ”".19A": J, 
mele 

5000 END IF 

S010 BEEP 

5020 PRINT 


5030 PRINT USING “{0X,"°"°NOTE: °".ZZ."" X-Y¥ oairs were stored in plot data file 


a 





Be 0a 3 J.P files 


5040 
3050 
5060 
S070 
5030 
SUE 


iF Inf={! THEN 

PRINT USING -16%.22, ° ° °,pairs of Nr-F are stored in file “",14A°;:J.Mrts 
END IF 

ASSIGN a@File TO ~* 

ASSIGN @Filel TO * 

ASSIGN @Filep TO * 

END 

DEF FNPvst(Tsteam) 

DIM K(8) 

DATA -7.691234564,-26.08023696,-168.1706546,64.23235504.-118.9646225 
NAIA) 4.167 11732,20.35750676,1.—E9.6 
READ K(*#*) 

T=(Tsteamt273.15)/647.3 

Sum=) 

FOR N=0 TO 4 
Sum=SumtkK (Nd *C 1-7) 7 CN 1) 

NEXT N 

Br=Sum/(T#( 14K (5) #C1-T)+K (6) *C1-T) 2) C1 -TI/ CK C7) C1 -T) UO 24K 08) ) 
Pr=EXP(Br) 

P=22120000Pr 

RETURN P 

FNEND 

DEF FNHfg iT) 

Hf a=2477200-2450#¢(T-10) 

RETURN Hfg 

FNEND 

DEF FNMuw(T) 

A=247 .3/¢(T+133.15) 

Mu=2.4E-5*107A 

RETURN Mu 

FNEND 

DEF FNVvst(Tt) 

P=FNPvst(Tt) 

T=Tt+273.15 

X=1500/T “a 
Fil=1/(t+T#1.E-4) 
F2=(1-EXP(-X))272.5*EXPCX)/X*.5 
B=.0015*F1-.000942F 2- .0004882*x 
K=2*P/(461.52*T) 
Y=C1+01+2*BeK)* .5)/K 

RETURN V 

FNEND 

DEF FNCpu tT) 
Cew=4.21120858-T#¢2.256826E-3-T*(4.42361E-5+2.71428E-/7«T)) 
RETURN Cpw*1000 

FNEND 

DEF FNRhow(T) 

Ro=999 .52946+T*¢« .01269-T#*(5.482513E-3-T*!.234147E-5)) 
RETURN Ro 

FNEND 

DEF FNPru(T) 


Meum Io7b6Uaa i Clee saocoe-2-1=(3.9608/5E-5-3.45025E-7*1)>) 


Pru=FNCpow( 7) *FNMuw( T)/FNKw¢ T) 
RETURN Pry 

FNEND 

DEF FNKu<T) 


mnw=noGe5G94+1* (2. 29645466-3-1*(1 .5097666-5-4.0581652E-8*T)> 


Pata anid) / 273.15 
Rut>. 922474K*(2 .8395-X¥¢ | .B007-X*¢ .52577-.07344*xX))») 
RETURN Ky 
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5640 
5650 
5660 
5670 
5680 
56390 
5700 
5710 
5720 
D7 oO 
5740 
5750 
5760 
5770! 
5780 
5730 
5300 
5810 


3840 
5850 
5860 
3870 
9880 
5890 
5900 
5910 
9920 
= Se) 
5940 
5950 
B00 
37 0 
9980 
5990 
6000 
6010 
6020 
6030 


SMEND 

DEF FNTanh(X) 

zeke (kD 

Q=EXP(-X) 

Tanh=(P+Q)/(P-Q) 

RETURN Tanh 

FNEND 

DEF FNTvsv(V) 

ebM 4Ce7 C7) 

Sum=C (0) 

FOR [=! TO 7? 
Sum=Sum+C (I) *V7I 

NEA la 

T=V*¢( .02635206956-V-=(9.7351313E-/7-Y<5.57680SE-11)> 
RETURN Sum 

FNEND 

DEF FNHf<T) 

Hf =T=(4.203849-1<«(5.88132E-4-1 «4.551603! 7E-6)) 
RETURN Hf*1000 

FNEND 

DEF FNGrad(T) 

Grad=37 .9853+.104388T 
RETURN Grad 

FNEND 

DEF FNTvsp(P) 

Tu=1!10 

i= 1-0 

Ta=(Tuttil)«.5 
Pe=FNPvst(Ta) 

IF ABS((P-Pc)/P)>.001 THEN 
IF Pc<P THEN Tl=Ta 

IF Pco>P THEN Tu=Ta 

GOTO 5910 

END IF 

RETURN Ta 

FNEND = 
DEF FNPvsv(V) 
P=8133-91393+2.23605!E+4=\) 
RETURN P 

FNEND 


a, 





PROGRAM TCAL 


iu ? 
eu)! 
(a 
130 
140 
150 
160 
170 
180 


meee NAME: 7 
Rev Sep: i) 


Coie Ce) C7 > 

Ditsemt< 10). 1109 20670) 

OomO UusbOgieeafa7 .94369,-75/7345.38259./73025595.81 
Dee -3247 46550 9,5.97668E11.-2.66192E12.3.94078E14 
READ C(~) 

CEA #09 


BEEP 

INPUT “ENTE MONTH, DATE AND TIME (MN:DD:HH:MM:SS)",BS 

J=0 

BUTPUT 70S: TD’ :BS 

OUTPUT 703. “TD” 

ENTER 709:AS 

oo USING "10X,.""Month. date and time = “",14A°%;AS 
Se 

INPUT “ENTER INPUT MODE ¢(1=30548, 2=FILE)” Im 

IF Im=1 THEN 

Boer 

INPUT "GIVE A NAME FOR DATA FILE",.D_fileS 

CREATE BDAT D_files.5 

Bo 

BEEP 

epee GIVE SNAME OF EXISTING FILE” D_files 


Beer 

INPUT “ENTER NUMBER OF DATA RUNS STORED" ,Nrun 
eNO IF 

BEEP 

(evie civE iA NAME FOR PLOT FILE’ .P_files 
CReatTEe BOAT P_ftiles.5 

ASSIGN sPlot TO P_files 

ASSIGN @File TO D_files 

IF Im=1 THEN 

Ber 


BEc 

INPUT “ENTER BATH TEMPERATURE”. T_bath 
OUTPUT 709:"“AR AFZO AL30” 

FOR I=0 TO 10 

OUTPUT 709:"AS SA” 

ENTER 709;:Emf C1) 

NEXT I 

OUTPUT 713:"“TIR2ZE” 

WAIT 2 

Spee 71337 1 

BUTPET 7133 °T2R2E" 

WAIT 2 

ENTER 713:1f2 

OUtreul Shite: T bath.emt ¢*),171.72 
ELSE 

ENTER '@Fiie:T_bath.&mfc*+). 11,72 
END IF 

J=J+1 

Dwa=0 

FOR I=0 TO 10 

TCL) =FNTvsv( ABS CEmf C1) )) 

D(C I)=T_bath-TCI) 

IF 1>4 THEN Dwa=Dwat)< I) 

NEXT I 

Dwa=Duwa/6 


CAL 
ecember 11, 1983 
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Dsa=(D(0)+DC1))=.5 
BDUTPUT @Plot:T_bath.Dsa. Dua 


PRINT 
PRINT 
PRINT 
PRU 
PRINT 
PRINT 
PRINT 
PRINT 
Bee 


IF Im= 


USING 
USING 
USING 
USING 
USING 
US ING 
USING 


1 TREN 


“10X.°° “Data set number = 


mix. Bath 1 (Cl) 


ei OG S0 ro. CO bath. |. 
"t0X." Thermocouple readinas (Deq CC): 
riers CoO 053k) on ee 1 C=) 
“10X.°°"Discrepancies (Deg C):""” 

PT Ae eiic -DD Tak ton 7 Gre 


BIG Ica) 


ee) ees 


(C) 


nd) 
14 


INPUT “ARE YOU TAKING MORE DATA (1=Y.0=N)?".Go_on 
pieces THEN 430 


IF J<Nrun THEN 430 


eNDae 


PRINT 


IF Im= 
PRINT USING “10X%.°"NOTE: 


SBS a 


PRINT USING “10X.""NOTE: Above analysis was performed from file "" 


Slo) 0S 
PrituleusiNG 5A. Plot data are stored in file 


ee 


ASSIGN #File TO « 
ASSIGN @#Plot TO «- 


END 


DEF FNIvsv(Emf > 
COMey Ge Ce C7>) 
Sum=C (0) 


FOR Is 


ites 


Sum=Sum+C C1) «ent “I 


NEX ie 


RETURN Sum 


FNEND 
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